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We generalize the static model by assigning a g-component weight on each vertex. We first choose a com- 
ponent (ff) among the q components at random and a pair of vertices is linked with a color n according to 
their weights of the component (/x) as in the static model. A (1 — /) fraction of the entire edges is connected 
following this way. The remaining fraction / is added with (q + l)-th color as in the static model but using the 
maximum weights among the q components each individual has. This model is motivated by social networks. 
It exhibits similar topological features to real social networks in that: (i) the degree distribution has a highly 
skewed form, (ii) the diameter is as small as and (iii) the assortativity coefficient r is as positive and large as 
those in real social networks with r reaching a maximum around / ~ 0.2. 

PACS numbers: 89.65.-s. 89.75.Hc, 89.75.Da 



I. INTRODUCTION 

Recently there have been considerable efforts to understand 
complex systems in terms of random graph, consisting of ver- 
tices and edges, where vertices (edges) represent individuals 
(acquaintances or their interactions) 01 IH EL Ell- In sucn com- 
plex networks, the emergence of a power-law degree distribu- 
tion, P(k) ~ k i, is an interesting feature. Such networks are 
called scale-free (SF) networks. To illustrate such SF behav- 
ior, many in silico models have been introduced, whose exam- 
ples include the Barabasi and Albert model 1 5], the Huberman 
and Adamic model |6[, etc. In those models, the number of 
vertices grows with time. 

The static model |7, 8] is another type of in silico model de- 
signed to generate SF networks, where the number of vertices 
is fixed. Each vertex is indexed by an integer i (i = 1, • • • , N) 
and assigned its own weight Wi — i~ a , where a is a tun- 
able parameter. Next, two different vertices (i, j) are selected 
with probabilities equal to normalized weights, Wi/ J2k w k 
and Wj/J2k w k> respectively, and are connected via an edge 
unless one exists already. This process is repeated until mN 
edges are present in the system, so that the mean degree is 2m. 
Then it follows that the degree distribution is SF with the ex- 
ponent 7 = 1 + 1 j a. Thus, tuning the parameter a in [0.5, 1), 
we can obtain a continuous spectrum of the exponent 7 in the 
range 2 < 7 < 3, for which the degree distribution has finite 
mean and diverging variance. Since the number of vertices 
does not grow, one may wonder if this model can be applied 
to evolving real world network. However, since the model 
network can be easily generated and exhibits little hump in 
the degree distribution, it can be useful to study many aspects 
of SF networks. 

In this paper, we generalize the static model by allow- 
ing a g-component weight (iu$ , , . . . , ) to each ver- 
tex i. We suppose that the /x-th component wf- of a ver- 
tex i represents its own weight or fitness to a subgroup (/i) 
(/i = 1 , . . . , q) in a society. For example, we suppose that 
two persons i and j are alumni of a high school, a subgroup 
(/i). They would have different weights w\ and Wj in the 
subgroup (/i), determined by their school activities. The per- 
son i and another person k are colleagues in a company, an- 



other subgroup [v). They have also different weights, 

and wj"\ by their positions in the company, the subgroup (y). 

Then the person i has weights and in different sub- 
groups, which are not the same in general. We make an edge 
between the pair in one color representing the subgroup 
(/i) and the pair (i,k) in another color for the subgroup (v). 
Vertices in the system are connected with edges in q different 
colors representing different subgroups. Subgroups are then 
connected each other by weak ties as explained later. Since 
our society comprises many different subgroups and a person 
can be acquainted with other people belonging to diverse sub- 
groups, this generalized static model is useful for modeling 
social networks. Meanwhile, it is noteworthy that our model 
is reminiscent of the generalization of the Ising model to the 
q-component cubic model [9] in equilibrium spin systems. 

So far, there have been many attempts to explain the struc- 
tures and the properties of social networks 1 10] . Recently, 
Watts et al. II ill introduced a hierarchical model for social 
network. In the model, individuals belong to groups that in 
turn belong to groups of groups and so on, creating a tree- 




FIG. 1 : A network of the g-component static model with parameters 
N — 80, m = 2, q = 4 and / = 0.2. Edges in four colors (red, 
yellow, green and blue) are the connections within each group. Edges 
in orange are those formed by maximum weights. Edges in more than 
two colors are colored in black. 
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like hierarchical structure of social organization. Here an in- 
dividual can belong to more than one group, as a result of 
which the distance between two persons is shorter than the 
ultrametric distance between them. Such hierarchical model 
illustrates well the small-world property of social network as 
implied in the Milgram's "six degrees of separation" fl^l . An- 
other simple social network model, introduced by Newman 
and Park Il3ll . is based on the concept of bipartite graph fl4ll 
and community structure 1 1 5 ] . This has an advantage toward 
explaining the non-trivial high clustering of real-world social 
networks. While our q-component model is similar to the hi- 
erarchical model and the community model in the spirit of 
dividing people into subgroups, however, we assign weights 
to each person for each subgroup, and connections are made 
following those weights. 

Social network exhibits an interesting feature in the degree- 
degree correlation function, different from biological or in- 
formation networks. Newman flrjll studied the degree-degree 
correlation in terms of the correlation function between the re- 
maining degrees of the two vertices on each side of an edge, 
where the remaining degree means the degree of that vertex 
minus one. He introduced the assortativity coefficient r, de- 
fined as 

r = \ ^2jK e jk ~ QjQk), (!) 

where ejk is the joint probability that the two vertices on each 
side of a randomly chosen link have j and k remaining de- 
grees, respectively, q^ is the normalized distribution of the 
remaining degree q% = (k + \)P(k + 1)/ ^ ■ jP{j)i an d 
<7g = J2k h 2< lk — Efc kqk} 2 - Interestingly, complex networks 
can be classified into three types, having r < 0, r « and 
r > 0, called the dissortative, the neutral, and the assortative 
network, respectively II 1611 . Most social networks are assor- 
tative as shown in TABLE I, while the Internet and the pro- 
tein interaction network are dissortative. While many in silico 
models have been introduced, most of them are neutral. Thus 
it would be interesting to introduce an in silico model having 
the assortativity coefficient as positive and large as empirical 
values, which would enable one to understand a basic mech- 
anism of social network formation. We will show that such 
assortative networks can be generated via the q-component 
static model. 



II. MODEL 

The g-component static model network is constructed as 
follows. Initially, TV vertices are present in the system, rep- 
resenting TV people in a society. Each vertex is assigned a q- 
component weight [vip , wf 1 , . . . , w^), where % is the ver- 
tex index. w\ ^' , the /i-th weight of a vertex i, is given as i\ ^ , 
where is the rank of the vertex i in the /i-th subgroup. We 
take {£i,fj,, . . . , (n,^} be a random permutation of the integers 
{1, . . . , N}. is also taken to be a real random number dis- 
tributed uniformly in the range [0.5, 1). In general, ranks of 
a person for different subgroups should be correlated in real 
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TABLE I: The size N, the mean degree (k), the diameter d, and the 
assortativity coefficient r for a number of social networks. 



society; however, we take them as independent in this work 
for simplicity. As the number of people N becomes large, 
the number of distinct subgroups q can increase in real world. 
Thus, we set q to be q = sN with s <C 1. Then 1/s is the 
average number of people belonging to one subgroup. 

Edges are connected as follows: First, among the q com- 
ponents, we choose a component /i at random. Second, 
we choose two vertices with probabilities equal to 

normalized weights, pf"' = / Y^k w k^ an( ^ Pj = 

Wj /Ylk w k ' an< ^ a tt acn an edge between them with the 
/i-th color unless an edge in that color exists already. Edge 
color is distinct for each component. Note that the pair 
can be connected via more than one edge in different colors. 
Edges in different component are distinguished by their own 
colors. The process of attaching such edges is repeated until 
(1 — f)mN edges are added to the system. / is a parameter 
between and 1. We will see that m is related to the aver- 
age degree. Since the component was chosen randomly, the 
number of edges in one color is (1 — f)mN/q on average. 

To construct a minimal model mimicking social relations, 
we need elements playing the role of "weak ties" lEOll . So, 
to take into account of social relationships among people 
having different backgrounds, we suppose that additional so- 
cial relationships are formed following the maximum weights 
among the q components each individual has. Let Wi >m = 

max(p ! - 1 ' ) , . . . ,pf^) be the normalized maximum weight of 
vertex i. Then two distinct vertices i and j are chosen with 
probabilities, Wi, m /J2k w k,m and Wj, m / J2k w k,m, respec- 
tively, and are linked by a new color different from the previ- 
ous q colors unless such an edge exists already. This process is 
repeated until fmN edges are formed. Edges formed by such 
maximum weights can be regarded as weak ties, introduced 
by Granovetter 1 20] which play an important role in social 
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TABLE II: Typical simulation results of the diameter d and the assor- 
tative coefficient r obtained under selected conditions of TV, m and 
q with / = 0.2. 
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FIG. 2: The degree distribution P(k) vs the degree k obtained with 
N = 10000 and m = 2 for q = 10, 100 and 1000. 
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FIG. 3: The assortativity coefficient r vs the parameter / for various 
m and q. Inset: / dependence of the diameter. 



networks, connecting different subgroups. We find that the 
assortativity coefficient is enhanced by the presence of such 
weak ties. 

Networks constructed in this way have mN edges with 
(g+1) colors representing internal structure of subgroups. So, 
some pair of vertices are linked by more than one edges in dif- 
ferent colors, albeit such incidences are not so frequent when 
q is large. However, when we measure the network proper- 
ties such as the shortest pathways, the degree of a vertex, the 
assortativity coefficient, and so on, we regard those multiple 
edges as a single one. Thus the mean degree (k) is slightly 
less than 2m by about 5% for typical networks we consider 
below. A small size network constructed in this way is shown 
in FIG. 1. 



III. SIMULATION RESULTS 

We perform numerical simulations for various values of q, 
/, m and A, and examine the diameter d and the assortativ- 
ity coefficient r as functions of those parameters. Here the 
diameter is the average distance between a connected pair of 
vertices along the shortest pathways. 

First, the shape of the degree distribution depends on the 
number of subgroups. For small q, the degree distribution fol- 
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FIG. 4: The assortativity coefficient r vs In In N for various m and 
q values. 



lows a power law with 7 w 2.8, however, for large q, it has a 
highly skewed form, approximately obeying a power law for 
a part of its range, and having an apparently exponential cut- 
off for larger k. The exponential cutoff for large k originates 
from the randomness of ranks of each subgroup, and the SF 
behavior for intermediate k does from the SF behavior of each 
subgroup. A similar crossover behavior can be found in real 
social networks, for example, the collaboration networks of 
physicists, biologists and movie stars l2~fll . 

Second, we examine the assortative coefficient r as a func- 
tion of / for a fixed A and several values of m and q. As 
shown in FIG. 3, the assortativity coefficient exhibits a peak 
around / w 0.2, meaning that the connections among sub- 
groups are mostly optimized. Thus we limit our further con- 
sideration to the case / = 0.2. Meanwhile, the diameter 
gradually decreases with increasing / as shown in the inset 
of FIG. 3. 

Third, we study the assortativity coefficient r as a function 
of A for various m and q. It is likely that r increases with 
increasing A as r ~ In In A as shown in Fig. 4, but it seems 
to saturate for larger A. It also increases with increasing m 
and q, as shown in FIG. 4. Thus the g-component static model 
exhibits r values as large as empirical values listed in TA- 
BLE I. Some numerical results of r listed in TABLE II show 
a quantitative agreement with the ones obtained in real social 
networks. 

Fourth, the diameter d is investigated as a function of the 
number of vertices A for various m and q, as shown in FIG. 5. 
It is found that the diameter is proportional to d ~ In A as in 
the case of random graph, in which d ~ lnA/ln(fc) l22ll . 
Furthermore, the diameter becomes smaller as m increases, 
which is like the case of random graphs. However, the di- 
ameter is almost insensitive to q. To test the so-called "six 
degrees of separation", we extrapolate the straight line in the 
semi-logarithmic plot of FIG. 5 to large A for m = 10 and 
q = A/100. We obtain d « 6.0 for A = 10 s and d w 6.7 for 
A = 10 9 , in reasonable agreement with the Milgram's "six 
degrees of separation" II 211 . Here the choice of m = 10 and 
q = N/ 100 comes from the facts that a person knows about 
20 other people on average (see Chapter 5 of Ref. 12311) , and 
there are about 100 members on average in a subgroup 1 24]. 
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FIG. 5: The diameter d vs the size N for the same parameters used 
in FIG. 4. 



the clustering coefficient behaves as C(N) ~ N 1 * 7 3 t)/(t 1 ) 
Q. Thus when 7 = 3, C(N) ~ AT" 1 . For the q- 
component static model, while r is not close to zero, the rule 
of attaching edges is such that no explicit degree-degree cor- 
relation enters, so that it is natural to expect the behavior 
C(N) ~ N^ 1 . Indeed, the measured behavior is close to 
the expected one as shown in FIG. 6, but the slope in FIG. 6 
exhibits small deviations for smaller q. For neutral networks, 
it is known that the clustering coefficient of a vertex with k 
is almost independent of degree k. Even in our case, we find 
that indeed C(k) is independent of k for different N as shown 
in the inset of FIG. 6. 



IV. CONCLUSIONS 
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FIG. 6: Plot of C(N) vs the size N for various m, q and / values. 
The slopes of the fit lines are —1.0 for solid and -0.9 for the dashed 
lines, drawn for the eye. Inset : Plot of C(k) vs degree k for various 
N with fixed m = 2, q = N/20 and / = 0.2. 



We have introduced the g-component static model assign- 
ing a g-component weight to each vertex. The weight of a 
given component of a vertex mimics a weight or fitness of 
that person in that subgroup. Through this model, we obtained 
the diameter of the acquaintance network as small as the Mil- 
gram's "six degrees of separation" and the assortativity coef- 
ficient as positive and large as empirical values for a variety of 
social networks. Moreover, we obtain the degree distribution 
in a skewed form, which is also similar to those of real world 
social networks. The clustering coefficients C(N) and C(fc) 
behave as those of a neutral network, being due to the absence 
of intrinsic degree-degree correlation. Such deficiency of the 
present model may be improved by introducing the hierarchi- 
cal structure among subgroups, or correlated ranks of a person 
for different subgroups. 
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Fifth, one of the properties well studied for complex net- 
works is the clustering coefficient C, which is defined as the 
average over all vertices of the ratio of the number of trian- 
gles connected to a given vertex to the number of triples cen- 
tered on that vertex. It is known that for the neutral network, 
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